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DISCLAIMER 


GPA publications necessarily address problems of a general nature and may be used by anyone 
desiring to do so. Every effort has been made by GPA to assure accuracy and reliability of the 
information contained in its publications. With respect to particular circumstances, local, state, and 
federal laws and regulations should be reviewed. It is not the intent of GPA to assume the dutics of 


employers, manufacturers, or suppliers to warn and properly train employees, or others exposed, 
concerning health and safety risks or precautions 


GPA makes no representation, warranty, or guarantee in connection with this publication and 
hereby expressly disclaims any liability or responsibility for loss or damage resulting from its use 


108s Or damage resi hg from ifs uw 


or for the violation of any federal, staic, or municipal regulation with which this publication may 


conflict, or for any infringement of letters of patent regarding apparatus, equi 
covered. 
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1 Scope 


1.1 This standard presents procedures 
for calculating, at base conditions from 
composition, the following properties of 
natural gas mixtures: gross heating value, 
relative density (real and ideal), and 
compressibility factor. The effect of 
water upon these calculations is 
somewhat complicated. Because this 
document relates primarily to custody 
transfer, the water effect included is an 
acceptable contractual calculation. The 
Appendix of this standard contains a 
detailed investigation of the effect of 
water and detailed derivations of the 
equations presented in the standard. 


2. Summary of Method 


2.1 From the composition (expressed. 
in mole fractions) of a natural gas 
sample, it is possible to calculate the 
gross heating value, relative density and 
compressibility (real gas) factor for the 
sampie. When analyzing the sample for 
composition, it is essential to include ail 
components (other than water) with moie 
fractions greater than or equal to 0.0001 
in the analysis (some routine analyses 
ignore constituents such as He and HS 
but ihey are important for accurate 
calculations). The gas sample should be 
collected according to the latest version of 
GPA Standard 2166 and the sample 


analysis should conform to the latest 


version of GPA Standard 2261 (or to 


other eatable methods) 
olicr acceptapie Meticds). 


3. Definitions 

3.1 Gross Heating Valuc - the amount 
of energy transferred as heat per mass or 
mole from the complete tdaat 
mole from the complete, idcal 


combustion of the gas with oxygen (from 


air), at a base temperature in which all 


water formed by the reaction condenses 
to hi 


uid (as exnlained in the Anpendix 
RC (as exprained In the Appendix 


this is a hypothetical state, but it is 
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acceptable for custody transfer). If the 
groce heating value has a volumetric 
rather than a mass or molar basis, a base 
pressure must be specified along with a 
base temperature. 


3.2 Relative density - the ratio of the 
mass density of the gas at the 
Measurement temperature and pressure to 
the density of dry air (the assumed 
composition of air appears in the 
Appendix) at the same temperature and 
pressure. Adjusting the result to the 
hypothetical ideal gas state converts the 
relative density into the molar mass ratio. 


3.3 Compressibility Factor - the ratio 
of the actual volume of a given mass of 
gas to its volume calculated from the 
ideal gas law using given conditions of 
temperature and pressure. 


3.4 Dimensions - 


1 Btu/lbm = 2.326 Jeg! (exact) 
1 Ibm = 453.59237 g (exact) 
1 Btu = 1055.0559 J 
1 atm = 101325 Pa (exact) 
= 14,69595 lbfein (psia) 


1 ft? = 0.0283168 m3 


4. Equations for Custody 
Transfer Calculations 


4.1] Gross Heating Value (Volumetric 
Basis) 


AY" (sat) = (t~ x, iv" (dry) (2) 

where Hv‘ is the gross heating value per 
volume al base temperature and pressure, 
superscript id denotes an ideal gas 
property, dry denotes dry gas, sat denotes 
gas saturated with water, x, are mole 
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fractions, N is the total number of 
components (excluding water) and x,, is 
the mole fraction of water in the gas. 
This Standard assumes that compositions 
are reported dry which is the usual case 
The Appendix illustrates procedures for 
the situation in which compositions 
include water. The latter quantity can be 
calculated from 


% = FTP, @) 


where P* is the vapor Pressure of water 
at the base temperature and P, 'b is the base 
pressure. The following table provides 
(1-x,) for some common base pressures 
used in the United States with a base 
temperature of 60°F where the vapor 
pressure of water is 0.25636 psia. 


B,(psia) 1-x, 

14,50 0.9823 
14.696 0.9826 
14.73 0.9826 
15.025 0.9829 


The next table presents (1-x,,) for some 
common base temperatures used outside 
the United States with a base pressure of 
1 atm. 


z,(°c} (1-x,) 
0 0.9940 
15 0.9832 
20 0.9769 
25 0.9687 


4.2 Relative Density 
G=dld,=(MPT.Z,)(M,P,TZ) (4) 
if T=T, and P= P, then 
G=(M/M,)\\Z,/Z)=G{Z,/Z) (5) 
where G is relative density, d is mass 


density, M is molar mass, P is pressure, 
T is temperature, Z is compressibility 
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factor and subscript a denotes a property 
of air, Calculation from Cumposition 
uses: 


GY =xG" +x,Gl+---+x,G! 


we ria 
KG, 


os 
oN 
a 


4,3 Compressibility Factor 


Al base conditions (near ambient) a 
simple expression which provides the 
compressibility (real gas) factor within 
experimental error for natural gas 
mixtures is: 


Z=1-# [San] (7) 


where the b; are “summation factors” as 


defined in the Appendix and listed in 
Table 1. Note: AGA-8-92 is the 
required compressibility factor 
method to be used in association with 
calculations under API Chapter 14,3 / 
ANSI 2530 / AGA-3 f GPA 8185. 
These two methods, in general, do not 
provide exactly the same value. The 
difference should be on-the-order-of 1 
part in 10,000 which is well within 
experimental error. 


De Example Calculations 


Si 
3.1 


G# and b, (along with heating values on 


mass and molar bases as well as net 
heating values) for the normal 


constituents of natural gas mutates 
required to caleuvlate the vai 


properties covered in this standard. Table 


2 and Table 3 provide calenlations 


illustrating custody transfer applications. 
The Appendix contains derivations and 


calculations illustrating the more nearly 
correct treatment of water. 
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6. Computer Program 


6.1 A listing of a Basic computer 
program which can perform eithcr 
custody transfer or engineering 
calculations is supplied with this 
standard, Typical output from this 
program precedes the listing. The values 
from computer output for some of the 
properties differ from those reported in 
the Tables (usually by one part in the last 
place). This is caused by different round 
off and truncation protocols between the 
computer and the hand calculator. 


7 Caution 


7.1 The properties reported in this 
document derive from experimental 
measurements which, in general, are 
accurate to no better than 1 part in 1000. 
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The extra digits. which appear in the 
Tables alleviale prublems assuciaied with 
round-off and internal consistency, but 
they are not significant. 


For numerical reasons, the values 
of some properties in this standard may 
differ slightly from those in GPA 
Standard 2145 (1 part in 10000 or Jess). 
These discrepancies are well within 
experimental error. In such cases, GPA 
2145 takes precedence for accounting 
purposes and GPA 2172 takes 
precedence for technical calculations. The 
values in GPA 2145 have been adjusted 
to achieve complete numerical 
consistency within that table. 
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APPENDIX 


DETAILS OF CALCULATION METHODS AND TREATMENT OF WATER 


Custody transfer of natural gas utilizes a simple pricing equation which states that the cost 
of gas is the rate of energy released upon combustion multiplied by the price of gas per energy unit 
multiplied by the time or accounting pericd. The rate of energy rcleased upon combustion is the 
product of the heating value of the gas and the flowrate of the gas. The flowrate of the gas requires 
knowledge of the compressibility factor and the relative density of the gas. All three custody 
uuansfer properties (heating value, compressibilily factur, and relative densily) can be calculated 
from the composition given pure component property tables. 


This appendix presents rigorous equations to calculate from composition the custody 
transfer properties of natural gas. The equations for calculating the properties of dry natural gas are 
well known, but this appendix also presents an account of the effects of water contained in the gas 
and in the air used to burn the gas. 


EQUATION DEVELOPMENT 


The heating value of a natural gas is the absolute value of its enthalpy of combustion in an 
ideal combustion reaction. The heating value is, therefore, an ideal gas property that can be 
calculated unambiguously from tables of pure component values and it has no pressure 
dependence. 


An ideal combustion reaction with fuel and air in the ideal gas state and the possibility of 
water in the fuel and air is: 


C,H yS,(id} + (a+ B14 + 7G + €)0, (id) + 0.04383 (a + B/4 + 7G + e)Arlid) 
+[0.00162(a 4+ B/44 y)(l4e)4 xo M1 —xy— x }ICO,(id) 
+ [3.72873(0:-+ 8/4 +y\(1 +e) + xy M1 xy — xc) IN (id) + (nf + 02 )H,O(id) 


= (a +0.00162(a + 8/44 y(t €)+ x, /(1- xy — x,)]CO,(id) 
+ A, H,Oid) +n), O(1) + SO, (id) 
+13.72873H(0c+ B/4 + y(t €)+ ay L(1— ay — x¢)IN, (id) 
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where "id" denoies the ideai gas state, a, 8 and y are stoichiometric coe! 


cients, £ is the fractinn 


= 


excess air, the composition of 


contained in the gas, 1) are the moles of water contained in the air, n%, are the moles of water 
contained in the product gas mixture, n,, are the moles of gas which actually condense, X¢ is the 


mole fraction of CO, in the gas and x, is the mole fraction of N, in the gas. If air has been 
injected into the gas, it is assumed that the effect is accounted for in the excess fraction ¢. Fuel gas 
mixtures would have non-integer values of a, 8 and y. 


di sed to “speciaior” water which enters the reaction carried by the 
gas or air). If we assume that the water formed in the reaction remains in the ideal gas state, the 
heat 8 termed “net.” if we assume that the water formed in the reaction condenses totally 
to the liquid state, the heating value is termed “gross.” The gross heating value is greater than the 


net heating value by the ideal enthalpy of vaporization for water; 


Heating value (gross) - Heating value (net) = H,,(id) -— H,,(1) (2) 


itis possibie to caicuiate a reai gas heating vaiue rather than using a hypothetical state, but 
the calculations are tedious, the numerical values are negligibly different and the mathematical 
implicity of the defining equation is fost. It is customary in the gas industry to use gross heating 


value for most calculations, so for the remainder of this appendix the term “heating value" refers to 
the gross value. 

Heating value is measured on a mass or molar basis and converted to the ideal gas state for 
reporting. Thus, at any given temperature the heating value is 


N 
An® = x.Hn/ Q) 
tal 


Hm" = iy x,M,Hm," J SM, (4) 


where Hn" is the heating vaiue in energy per moie, x, is the mole fraction, N is the number of 
components in the mixture, Hm is the heating value in energy per mass, and M is the molar 


mass. Clearly, Hm" multiplied by the molar mass (with units of mass per mole) gives Hn”. 


Both Hn and Hm are independent of pressure, but both are functions of temperature and 
composition. 


10 


Copyright por American Petroleum Institute 
Sat Nov 03 22:10:45 2001 


GPA STD*2172 86 MM 3624699 OO1LS47? 5b? 


The natural gas industry uses heating value with dimensions of energy per volume in its calculations. 


These dimensions result from multiplying Hin or He# by density or mass density of the Ideal gas 


Hr, ,i¢ 


mayo= 


i=l 6) 
where Av" js the heating value in energy per volume, P is absolute pressure, T is absolute temperature and R 


is the gas constant (=8.314471 Jemol-!*K-1=10.7316 psia*ft3*Ibmol-t*"R-!). Hv" depends upon temperature, 


. Table A2 contai 


s values for Hv“ at 60 °F and 14.696 nsia. These valnes 


valid at the specifie 


the new Pand lA 
te new F and 14, 


Hy" = Hy"(Table2,P/14.696] 6) 


When using Equation 6, Hv" (Table A2 ) should be calculated using the values from Table A2 in Equation 5; the 
individual values in Table A2 should not be converted. Conversion to another temperature is more 
complicated. Heating value data exist at 25 “C based upon the reaction: 


C, HS, (id) + (0+ B14 + ¥)O, (id) = CO,(id) + (Bi DH D + BO, Cd) n 


The experiments use pure oxygen and are corrected to stoichiometric proportions. It is necessary to correct the 
sensible heat effects to arrive at a different temperature: 


7 


r 71 
Hn'*(T)= Hnt(Q5)+ |). |= ca aes PE 


7 7 (8) 
where 
id ytd 4 il 
LG = AC pico, +B! Ho FIGs % 
id _ pie Ace 
Oe = Chains, $+ B14 + PICK, as 
Ge. , : 
and ~“? is the ideal specific heat at constant pressure, r denotes reactants and rr denotes products. 
The cost of gas comes from the simple accounting equation 
c=O%p"At (1) 


11 
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where ¢ is the cost, Q" is the ideal rate of energy transfer, p“ is the price of gas per ideal energy 


unit and Az is the accounting period. Using real gas rate of energy transfer merely req e 
of gas per real energy unit which would differ from that in Equation 11 in exact proportion to the 


O"p" = Op (12) 
Q" results from multiplication of heating value by gas flowrate 


QO" =AHn" = mHm® = V4 Ay" (13) 


al ely. Gas industry 
practice dictates use of real gas volumetric flowrate most t flowmeters, such as orifices, provi le 


naturally the mass flowrate which, if used, would lead to more accurate me: 


necessary to convert the real gas flowrate into an ideal gas flowrate to use in Equation 13 


id . 
where #, rie and V“ are the molar, mass and ideal g 


V=v¥z (14) 


where Z is the compressibility factor (which is defined as the ratio of real gas volume to ideal gas 


volume). Now the energy flowrate becomes 
OY =(V/ Zz) Hv (15) 


The factor //Z in Equation 15 tigorously converts the real gas flowrate into an ideal gas flowrate. 


Tt does not convert heating value into a real gas property. Often calorimeter and chromatograph 
manufacturers report the value of Hv /Z as output. This is a convenience for the user allowing 
inmicdiate multiplication by V and thus satisfying Equation 15. 


‘The truncated virial equation of state satisfactorily represents Z at pressures near ambient 
Z=14+ BPI RT (16) 


where B is the second virial coefficient which is a function only of temperature and composition. 
The virial coefficient is (to a good approximation) 


where the 4, are summation factors" which equal 


i2 
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(18) 


Table 2 lists values for b,. 


Another property required to evaluate flowrate is the molar mass of the gas. Historically, 
the gas industry obtains this value from measurements of the gas relative density which is the 
mass density of gas divided by the mass density of air 


G=d/d,=MPT.Z,/M,PTZ (19) 


where dis mass density and subscript « refers to air. If the P and T of gas and air are identical (as 
recommended for measurement} 


G=(MIM,\Z,/Z)= G"Z,/Z (20) 


where G” is ideal relative density which equals the molar mass ratio of gas to air. The molar 


mass of air for the assumed composition is 28.9625 g*mol-!, G4 is a simple function of 
composition 


N 
GY =} x.G" @Q1) 


Table A2 lists values for G,". 


13 
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ACCOUNTING FOR WATER 


If the gas contains water (or must be assumed to be saturated) but the compositional 
analysis is on a dry basis, it is necessary to adjust the mole fractions to account for the fact that 
water has displaced some gas, thus lowering the heating value. The mole fraction of water in the 
gas results from the definition of relative humidity 


= Apt po - “y 
x, = WPS IP =n, (l+n,) (22) 
(based upon one mole of the fuel C,27,5,}) whore 


vapor pressure of water and x, denotes moles of water. For saturated gas h* is unity. 
Rearranging Equation 22 gives the moles of water 


n, =x, /(1-x,) (23) 
The corrected mole fractions then become 
1 | 1 


ep - | legicy. x; (24) 
[itn | ‘Lies, x] 


and the heating value becomes 


WN 
Hy =(i = > xi" Ay (25) 


isl 


where water is not included in the N components of the summation. If the compositional analysis 


determines x, and water is included in the WN components of the summation 


gy id 
ne 


a condensation effect. Water carried by wet gas (spectator water) dees not aeuially con 


by definition, water formed in the reaction contributes to heating value. 


Accounting for water in the above manner is sufficient for defined custody transfer 


s, but when trying to model actual situations the question becomes much morc 


complicated. It is obvious that all of the reaction water actually cannot condense because in a 


7 tr are tom wirntas eotuestes thy odu: 
situation in which both gas and air are dry some of the reaction water saturates the product gascs 


14 
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and the remainder condenses. It is possible to account for these effects in a general manner. To do 
So, it is necessary to calculate n’, ni, nm’, and nj. 
ne Ait (xy 4 xe)i(i— Ay — Xp) tne] = hEPo IP 
ng = WPS PAG — ay — xg i MPS AP) (27) 
ne “MATTIAS + Bi4+y\0 tejtnij=h*PrTrP 
ni = 4.774 18(er + 8/4 + y\I + eh’ POs P\(1— “Po / P| (28) 
WHat y + (xy +xc)1= xy — Xe) + (@+ B/4+ [0.001 62(1 + £) 
+3.72873(8 + €) + 0.04383(1 +2) +e]+n'} = PosP 
Baty t(xy +x,)A1- xy — x) + (@ + B/4 + y)[0.00162(1 + €) 
+3.72873(1 + €) + 0.04383(I + €) + €]}(P? (P)(1- P27 P) (29) 
nm, = B/2+n8 +n nt, G0) 


where 4, is the relative humidity of the air. Equation 27 and 28 are reformulations of Equation 22 
to reflect inlet conditions. Equation 29 reflects Equation 22 for the saturated product gas (it must 
be saturated before any water can condense). Equation 30 is a water balance: $ /2 are the moles 


of water formed by the reaction, n? +n’, are the moies of water which enter with the gas and air, 


m, are the moles of water which saturate the product gas and a! are the moles of water which 
condense. Therefore, the complete correction for the effect of water on heating value is 


Hy" = Hy" (Equation25or26) + {(h*P° | PY — xy —x¢)(1—ASPT IP) 
44.77418(01 + B/4 +) + eh Po / PY — hh? / PY 


Tal RtAG 


ARTY + hy — XC) — dy — a) 4 (G+ B/44 Y\B.TIAIB 4 4.7741 88)] 
p 


Depending upon the relative humidities of the gas and air, the observed heating value can be 
greater or smaller than that calculated using Equation 25 or 26. A humidity of air exists for each 
gas above which Mv“ is greater than that calculated by Equation 25 or 26. That critical value 


15 
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ain a Babe dat oescn, oats Ber 
the amount of excess air. For 


pure, dry methane with no excess air h, = 0.79345. 


REAL GAS PROPERTIES 
In principal, we have enough information to convert the heating value to a real gas property 


{it is not necessary to do so for relative density because the molar mass ratio, G“, is the desired 
property). This is simply a matter of evaluating the integral: 


: ra @2) 


where 
(Sf) =v-7[ 2) =5-1@ rer (33) 
AP oh \ GE Jp ai Gi 


oe 


T ine temperature dependence of 2 must be defined, but in the 
custody transfer Tegion it is easy to do so, The products and reactants again correspond to 


where) 
whereY is the molar vo 


po Pee oar 


itis obviously no to make the required calculations to convert the heating value 
into a real as property, it serves n custody transfer Purpose 1 to do so. As we have seen, the cost 
equation is unchanged; the calculations, while obvious, are tedious; Hn is slightly different from 
Hnid pecale th the base pressure is low; the likelihood of having all the information required to 


joy Equation 31 is remoie. 


@ 


The heating value is defined in a hypothetical state. it is not possibie, at base conditions, to 
have all the water formed in the reaction be either all Ba or all liquid; some of the water formed is 
Thies 


in cach state. Thus, if the definition is of a hypothetical state, using a hypothetical real gas rather 
than an ideal gas state adds nothing but complexity. 


16 
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Table Al 
Assumed Composition of Air from: Jones, F. E.; J. Res. Nat. Bur. Stand.; 83; 419 (1978) 
Component Mole Fraction 

No 0.78102 

On 0.20946 

Ar 0.00918 

COz 0.00034 

1.00000 
Table A2 
Properiies of Naiurai Gas Componenis at 60 °F and 14.696 psia 
ion Me Dorada git Bb (psig?) Hy (Brae f°) 
Component rormua aod Ne e be aed 
Water H,0 0.62202 0.0623 50.312 
Nitrogen Ng 0.96723 6.0044 0 
Carbon Dioxide CO» 1.5196 0.0197 0 
Air Table 1 1,00000 0.0050 0 
Methane CHy 0.55392 Q.0116 1910.0 
Ethane CoH 1.0382 0.0239 1769.7 
Propane C3Hg 1.5226 0.0344 2516.1 
iButane Cato 2.0068 9.0458 3251.9 
Butane C4Hi¢ 2.0068 0.0478 3262.3 
i-Pentane CsH12 2.4912 0.0581 4000.9 
Pentane CsHy2 2.4912 6.0631 4008.5 
Hexane CéH14 2.9755 0.0802 4755.9 
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NATURAL GAS PROPERTIES FROM COMPOSITION 


CONDITIONS: T(heating value) = 60 DEG F T(density) = 60 DEG F P= 14.696 psia 
ANALYSIS 

x(CH4) = 0.8302 x(C6l114) = 0.0030 x(H20) = 0.0000 
x(C2H6) = 0.0745 x(C7H16) = 0.0000 x(CO) = 0.0000 
x(C3H8) = 0.0439 x(C86H18) = 0.0000 x(N2) = 0.0032 
x(iC4H10) = 0.0083 x(C9H20) = 0.0000 x(02) = 0.0000 
x(C4H10) = 0.0108 x(C10H22) = 0.0000 x(H2S) = 0.0000 
xGC5SH12 = 0.0031 x(H2) = 0.0000 x(Ar) = 0.0000 
x(C5SH12) = 0.0025 x(He) = 0.0003 x(CO2) = 0.0202 
Molar Mass = 20.248 GROSS HEATING VALUES 
Molar Mass Ratio = 0.69909 22111 Btuelb! 

Relative Density = 0.70109 1179.7 Btus(Ideal CFY! 
Compressibility Factor = 0.99677 1183.6 Btue(Reat CF)! 


SATURATED ANALYSIS 


x(CHA4) = 0.8157 x(C6H14) = 0.0029 x(H20) = 0.0174 
x(C2H6) = 0.0732 x(C7H16) = 0.0000 x(CO) = 0.0000 
x(C3H8) = 0.0431 x(C86H18) = 0.0000 x(N2) = 0.0031 
x(iC4H10) = 0,0082 x(C9H20) = 0.0000 x(02) = 0.0000 
x(C4H10) = 0.0106 x(C10H22) = 0.0000 x(H2S) = 0.0000 
xGiC5H12 = 0.0030 x(H2) = 0.0000 x(Ar) = 0.0000 
x(CSH12) = 0.0025 x(He) = 0.0003 x(CO2) = 0.0198 
Molar Mass = 20.209 GROSS HEATING VALUES 
Molar Mass Ratio = 0.69775 21767) Bulb"! 

Relative Density = 0.70001 1159.2 Btue(Ideal CFy! 
Compressibility Factor = 0.99639 1163.3 Btu*(Real CF)! 
WET ANALYSIS 

CONDITIONS: 

RELATIVE HUMIDITY=0.5 HUMIDITY OF GAS=0.3 FRACTION OF EXCESS AIR=0.2 

SATURATED ANALYSIS 

x(CH4) = 0.8259 x(C6H14) = 0.0030 x(H20) = 0.0052 
x(C2H6) = 0.0741 x(C7H16) = 0.0000 x(CO) = 0.0000 
x(C3H8) = 0.0437 x(C86H18) = 0.0000 x(N2) = 0.0032 
xGC4H10) = 0.0083 x(C9H20) = 0.0000 x(O2) = 0.0000 
x(C4H10) = 0.0107 x(C10H22) = 0.0000 x(H28) = 0.0000 
xGiC5H12 = 0.0031 x(H2) = 0.0000 x(Ar) = 0.0000 
x(CSH12) = 0.0025 x(He) = 0.0003 x(CO2) = 0.0201 
Molar Mass = 20.235 GROSS HEATING VALUES 
Molar Mass Ratio = 0.69868 21919 = Btuelb1 
Relative Density = 0.70075 1168.8 Btue(Ideal CF)! 
Compressibility Factor = 0.99666 19 1172.7 Btu¢(Real CF)! 
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REM tA AR AA RARER EEE DECLARE fe RH RR RAR ae 


REM ft **#¢4 a Ate R ER H YARTABLES 6 RAR RR RK EER IR IKKE KIRKE ERE 


REM 
REM Al, AZ, A3 = SPECIFIC HEAT CONSTANTS 
REM AA, AB, AC, AD = CALCULATION VARIABLES FOR SPECIFIC HEATS 
REM ALPHA = NUMBER OF CARBON ATOMS IN THE FUEL 
REM B = SUMMATION FACTOR 
REM BETA = NUMRER OF HYDROGEN ATOMS IN THE FUEL 
REM BP = SECOND VIRAL COEFFICIENT 
REM Cow = CHECK VARIABLE FOR MOLE FRACTIONS SUM 
REM =A = FRACTION OF EXCESS AIR 
REM FS, FILENAMES = CHARACTER STRING VARIABLES 
REM FAC = CALCULATION VARIABLE FOR Z% 
REM GAMMA = NUMBER OF SULFUR ATOMS IN FUEL 
REM GI = MOLAR MASS RATIO 
REM GR = RELATIVE DENSITY 
REM Hl - HS, HC = CALCULATION VARIABLES FOR HEATING VALUE 
REM H = HEATING VALUE AT 25 C KJ/MOL 
REM HA = RELATIVE HUMIDI 
REM HG = RELATIVE HUMIDITY OF GAS 
REM HM, HV, HVR = GROSS HEATING VALUE 
REM HN = HEATING VALUE OF GAS 
REM HNA = DUMMY VARIABLE FOR HN 
REM ITR, ITR2 = COUNTER VARIABLES FOR LOOPS 
REM NAS = NAME STRING 
REM NAW = MOLES OF WATER ENTERING WITH AIR 
REM NGW = MOLES OF WATER ENTERING WITH GAS 
REM NVW = MOLES OF WATER IN VAPOR PHASE 
REM P = PRESSURE 
REM PC = CRITI P E 
REM PP = PRESSURE IN PASCALS 
REM Pw = VAPOR PRESSURE OF WATER 
REM RM = MOLAR MASS 
REM SUMK = SUM OF MOLE FRATIONS 
REM TL = BASE TEMPERATURE FOR RELATIVE DENSITY 
REM T2 = BASE TEMPERATURE FOR HEATING VALUE 
REM TA, ‘IB = DUMMY VARTABLES FOR T1 AND T2 
REM TC = CRITICAL TEMPERATURE 
REM TR = REDUCED TEMPERATURE 
REM VC = CRITICAL VOLUME 
REM xX = MOLE FRACTION 
REM XW = WATER SATURATED MOLE FRACTIONS 
REM Z = COMPRESSIBILTY FACTOR OF GAS 
REM 21, 22 = CALCULATION VARIABLES FOR WATER VAPOR PRESSURE 
REM ZA = COMPRESSIBILTY FACTOR OF AIR 
REM 
REY COCO GGG ORE OOOO ISITE ICCC GEIS SIE IEEE OIE IE III IIIS I I ie ie 
REM 
REM 
REM 
20 
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RR RR SEMIN UP BILE PAGE OO eI mer mar I WIC IE IH 


CLS : PRINT : PRINT : PRINT 
PRINT : PRINT SPC({20); "CALCULATION OF GROSS HEATING VALUE," 


PRINT : PRINT SPC(16); "RELATIVE DENSITY and COMPRESSIBILITY FACTOR" 
PRINT : PRINT SPC(26}); “for NATURAL GAS MIXTURES" 
PRINT : PRINT SPC(25); “from COMPOSITIONAL ANALYSIS" 


PRINT : PRINT SPC(36); "1993" 
PRINT : PRINT SPC(25}); "GAS PROCESSORS ASSOCIATION" 
PRINT : PRINT SPC(27); "6526 EAST 60th STREET" 


PRINT SPC(2S}; "TULSA, OKLAHOMA 74145 ¥. 


PRINT ULSA, OKLAHOMA 74145 USA’ 


PRINT : PRINT SPC(2S)}; 


PRINT : PRINT SPC(31}; “developed at" 

PRINT SPC(10); “Texas A & M University, Collega Station, Texas 77843 USA" 
PRINT : PRINT : PRINT : PRINT SPC(5); "{press the space bar)" 

10 I$ = INKEY$: IF Ig = "" GoTo 10 

CLS : PRINT : PRINT : PRINT : PRINT : PRINT 

INPUT "DO YOU WISH TO SAVE THE OUTPUT TO A FILE"; FS 

IF FS = “Y¥" THEN 20 

IF F$ = "y" THEN 20 

GOTO 30 

20 PRINT "Name the drive and file." 


PRINT "For example, B: STREAM] .DAT" 


INPUT "What are the drive and filename”; FileName$ 

OPEN "o", 1, FileNames 

REM 

REM *#*# 84 tee DIMENSIONING VARIABLES 08 00 IOGEAR ERR RH 
REM 

30 DIM X(21), TC(21), VC(21), A1(22), A2(22), R322) 

DIM RM({21}, NA$(21), B(21), H(21) 

DIM XA{21) 

REM 

REM *4*4 44 Rew WH READING IN DATA 444 200K KAKA AKAD ARR AIR IA EERE RRR RE 


REM 
FOR A =1 TO 21 

READ NAS {A) 
NEXT A 
DATA "(CH4) ", "{C2H6) ", "(C3H8) ", "(204H10) ", "(C4H10) ", "(iC5H12) 
DATA "(C5H12) =", "({C6H14) ", "(C7H16) “, “(C8H18) ", "(C9H20) “, “(CLOH22) 
DATA "(H2) “, “(He) ", *(H20} ", "{CO} ", "(N2) “, "(02) 
DATA "(H2S} ", "(AR) ", * (C02) * 


FOR A = 1 To 21 
READ TC{A}), VC(A), AL(A), A2(A), A3(A), H(A}, RM(A) 


NEXT A 

READ Al{22}, A2{22}, A3(22} 

DATA 190.55, .0000990, 4.1947, .003639, .00001490, 890.630, 16.0430 
DATA 305.33, .00014706, 5.9569, ,013770, 00001690, 1560,690 20,0700 


DATA 369.85, .0002000, 8.2671, .022860, .00001900, 2219.170, 44.0970 
DATA 407.85, .0002590, 10.8240, .031530, .00000820, 2868.200. 58.1230 
DATA 425.16, .0002550, 11.1090, .028750, .00001820, 2877.400, 58.1230 
DATA 460.43, .0003060, 13.4120, .035400, .00001400, 3528.830, 72.1500 
DATA 469.70, .0003040, 13.5870, .032880, .00001400, 3535.770, 72.1500 
DATA 507.40, .0003700, 16.1340, .039860, .00003600, 4194.950, 86.1770 


21 
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DATA 540.20, .0004320, 18.6420, .O4786C, 60003840, 4853.430, 160.2040 
DATA 568.83, .0004920, 21.1920, .054800, 00004300, 5511.800, 114.2310 
DATA 594.64, .0008480, 23.7300, .061720, .00004760, 6171.150, 128.2590 
LATA 617.60, .U0d6030, 26.280, .UoeYUY, .bUUUSZSU, GbZy./U0, 142.2850 
DATA 33.20, .0000650, 3.4330, .001550, -.00000740, 285.830 2.0159 
DATA 5.19, .0000573, 2.5000, .000000, .00000000, 0000.9000, 4.0026 
DATA 647.14, .0000559, 4.0260, .000410, .00000255, 44.106, 18.0153 
DATA 132.912, .0000930, 3.5030, .000090, .00000090, 282.980, 28.0100 
DATA 126.20, .0000891, 3.5020, .000060, .00000000, 0000.000, 28.0134 
DATA 154.58, .0000733, 3.5200, .000440, .00000280, 0000.000, 31.9986 
DATA 373.40, .0000982, 4.0700, .001180, .00000280, 562,010, 34.0800 
DATA 150.86, .0000746, 2.5000, .000000, .00000000, 0000.000, 39,9480 
DATA 304.21, .0000942, 4.3240, .005800, -.00000650, 0000.000, 44.0100 
DATA 4.797, 9043900, 9.0000013 

KS = "Y¥": LS = "Y" 

REM 

REM ****** INPUTTING TEMPERATURE, PRESSURE, AND MOLE FRATIONS *****s#e4 ket e oe 
REM 

40 CLS ; PRINT : PRINT PRINT ; IF K$ = "Y” THEN 50 

IF K¢ = “y" THEN 50 

PRINT : PRINT : PRINT IF KS <> "¥" THEN 80 

50 PRINT SPC(3}; "Do you use the same base temperature for density and heating 
value" 

INPUT Y¥$: IF Y$ = “Y¥" THEN 60 

IF y$ = ‘y" THEN 60 

INPUT "What is the base temperature for density? ", T1: TA = TL 

INPUT "What is the base temperature for heating value? ", T2: TB = T2 
GOTO 70 

60 INPUT "What is the base temperature? ", Ti: T2 = Tl 

TA = Ti: TB = T2 

70 INPUT "What is the base pressure? ", P 

Tl = TA: T2 = TB 

ITR2 = 0 

IF L$ = "¥" THEN 90 

IF LS = "y" THEN 90 

IF L$ <> "Y" THEN 100 

is LS <> "y" THEN 100 

90 FOR A= 1 TO 21 

PRINT "X"; NAS(A); "=": 

INPUT X(A) 

XA(A) = X(A) 

NEXT A 

100 FOR I = 1 TO 21 

K(I) = KA(I) 

NEAT I 

CLS PRINT : PRINT PRINT 

PRINT SPC(21}; “NATURAL GAS PROPERTIES FROM COMPOSITION": PRINT : 

TF Tl >= 30 THEN 110 

PRINT "CONDITIONS:"; SPC(1}; "Tlheating value} -"; T2; "C;"; SPC(2}; "Pldensit 
="; T1; "C;"; SPC(2); "B PB; 

T1 = T1 + 273,15: T2 = 12 + 273.15 

PP = P * 1000 

GOTO 120 

110 PRINT "CONDITIONS:"; SPC({1); "T(heating value) pened “3; SPC(2); 
"T({density) ="; Tl; "F; SPC(2); “P ="; P; “PSIA*: PRINT 

VTi = (TL + 459.67) / 1.8: T2 = (P2 + 459.67) / 1.8 

PP = P f .0001450377# 

120 PRINT SPC(35}; "ANALYSTS": PRINT 
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REM 


REM *+#s4¥eeeeseeaesees CALCULATING SUMMATION FACTORS ttt thteesaaa ea eanenne 
REM 

PUK 1 = 1] wu al 

TR = TC(I) / TL 

BP = VC(I) * (2.058978# * TR - 7.020596# * TR * 2 + 5.953652# * TR * 3 - 
2.057299% * TR * 4) 

B(I} = SQR(ABS(BP) / (8.314471 * T1)) 

NEXT I 

REM 


REM }*o 44 eee ee ee dees eee CALCULATING ZA Serb t ER RR ah hk te 


REM *44*44¢4 edd HE EX CEYTING ITERATION CHECK * AK KKK HERR ER RRA RARER ERA R EAR 


REM **#44¥44e2eeoee+ee* CHECKING MOLE FRACTION SUM *4*** te te RR Ei 


1 TO 21: SUMx = SUMX + X(I): NEXT I 

COW = ABS(SUMK ~- 1) 

IF CoW < .000005 THEN 130 

CLS : PRINT "WARNING: mole fractions da not sum to unity." 
PRINT : PRINT "Do you wish to re-enter your data?" 

INPUT YS 

IF Y$ = "Y* THEN 50 

IF YS = “y" THEN 50 

130 IF F§ = "Y" THEN 140 

If FS = "y" THEN 140 

GOTO 1760 

140 PRINT #1, SPC(21); "NATURAL GAS PROPERTIES FROM COMPOSITION": PRINT #1, 


PRINT #1, 


IF TA > 30 THEN 150 


DO 


PRINT #1, "CONDITIONS:"; SPC(1); "OQ;"; SPO(2); 
"T(density) ="; TA; "C;"; SPC(2); "PB P; "KPA": PRINT #1, 

GOTO 160 

150 PRINT #1, "CONDITIONS:"; SPC(1); “T(heating value) ="; TB; "F;"; SPC(2); 
"[Didensity) ="; TA; "F;"; SPC(2); "P =*; P; "PSIA": PRINT #1, 

160 PRINT #1, SPC(35); "ANALYSTS"; PRINT #1, 

REM 

REM ***4* CALCULATING MOLAR MASS, RELATIVE DENSITY, AND HEATING VALUE ***** 
REM 

170 GI = 9: HN = 0: RM = 90 

FOR I = 1 70 21 

TF X(T) = 0 THEN 270 

IF I> 4 THEN 180 

AA = I * AL(21) + (1 + 1) * A115) - ((3 * 1 +1) / 2) * at(18) 

AB = I * A2(21) + (1 + 1) * A2(15) - ((3 * 1 +1) / 2) * A2(18) 

AC = I * A3B(21) + (1 + 1) * A3(15) - ((3 * 1 +1) / 2) * A3(18) 
AD=1I+1 

GOTO 260 

180 IF I <> 5 THEN 190 

AR = 4 * Al(21) +5 * AL(15) - 6.5 * A1(18) 

AB = 4 * A2{21}) +5 * A2(15} - 6.5 + A2\18)} 
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AC = 4 * A3(2Z1) +5 * ASi1b) - 6.5 * A3(18) 

AD = 5 

GOTO 260 

190 IF 1 < 6 OR L > 7 THEN 200 

AA = 5 * A1(21) + & * A1(15) g * BA1(18} 

AB = 5 * A2(2Z1) + 6 * A2(15) - 8 * A2(18} 

AC = 5 * A3(21) + 6 * A3(15) - 8 * A3(18} 

AD = 6 

GOTO 260 

200 IF I < 8 OR I > 12 THEN 210 

AA = (I - 2) * Al(21} + (I - 1) * al(is) - (3 * I - 5) 7 2) * AL(ig) 
AB= (I - 2) * A2(21) + (1 - 1) * A2(15) - ((3 * I - 5) / 2) * AZ(i8) 
AC = (I - 2) * A3(21) + (1 - 1) * a3(15) - (43 * 2-5) 7 2) * A3(18) 
AD=I-1 

GOTO 260 

210 IF I <> 13 THEN 226 

ABR = A1(15) - .5 * Al(18) 

AB = A2(15) - .5 * A2(18) 

AC = A3(15) - .5 * A3(18) 

AD=1 

GOTO 260 

220 IF I <> 15 THEN 230 

AA = A1(15) 

AB = A2(15} 

AC = A3(15) 

AD=1 

GOTO 260 

230 IF I <> 16 THEN 240 

AA = Al(21) - .5 * Al(18) 

AB = A2(21) - .5 * A2(18) 

AC = A3(21) - .5 * AB(18) 

AD = 0 

GOTO 260 

240 IF I <> 19 THEN 250 

AA = A1(15) + A1(22) - 1.5 * aAl(13) 

AB = A2(15}) + A2(22) - 1.5 * A2(18} 

AC = A3(15) + A3(22) - 1.5 * A3(18) 

AD=1 

GOTO 260 

250 IF I <> 14 AND I <> 17 AND I <> 18 AND I <> 20 AND T <> 21 THEN 260 
AA = Al{I) 

AB = A2(I) 

AC = A3(I) 

AD = 0 

260 H3 = T2 - 273.15 

H2 = AD * (-.04135 * (H3 - 25) + 3.724E-05 * (H3 - 25) * 2) 
He = (AC AZ(Z}) * (H3 * HB + H3 * 50 4+ 625) / 3 

H4 = AA - Al(T) + (AB - A2(I)) * (H3 + 25) / 2 - HS 


HL = (3 - 
HC = H(I} - 


} 
25) * H4 * 8.314480000000001D-03 
Hl + H2 


AN = HN + X(I) * HE 
AN = HN - X(15} * H(15) 


HNA = HN 
270 NEXT I 
280 FAC = 0: 


RM = 0 


FOR I = 1 TO 21 


FAC = FAC + 
M 


RM = RM + X 


x(T) 
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TI 
Z=i1- PP * FAC * FAC 
HM = 1000 * HN / RY 


RM / 24.9025: GR = GL 


REM tebe EER RRR EEE PRINTING RESULTS (CCGG IO IER I Ix at 
REM 

1 FORI=17T07 

PRINT SPC(7); USING "x\ \= BREE x\ \= #.#8e# x\ \= 

4. eHee"; NAS(I); K(I); NAS(I + 7); K(I + 7}; NAS(I + 14); X(T + 14) 

NEXT I 

PRINT : PRINT SPC(7); USING “Molar Mass = FTF GROSS 
HEATING VALUES "; RM 
IF FS = "¥" THEN 290: IF FS = "y" THEN 290 

GOTO 45UU 

290 FOR I = 1 TO 7 

PRINT #1, SPC(7); USING "x\ \s 8. oe88 x\ \= #.ddd8 x\ 
\= B.#eRE "; NAS(I); X(I); NAS(I + 7); X(I + 7); NAS{I + 14}; X(I + 14) 

NEXT TI 

300 IF TA >= 30 THEN 320: IF TB >= 30 THEN 320 

PRINT SPC(7); USING "Molar Mass Ratio = FLARE fH 
k7/kg”; GI; HM 

PRINT SPC(7); USING “Relative Density = tHE E HE 
kJ/Ideal CM"; GR; HV 
PRINT SPC(7)}; USING "Comp 
kJ/Real CM"; Z2; HVR 

IF F$ = "¥" THEN 310 

IF F$ = "y" THEN 310 
GOTO 350 

310 PRINT #1, : PRINT #1, SPC(7); USING "Molar Mass = fH. 
GROSS HEATING VALUES"; RM 


ressi ctor = Bete? WHE 


PRINT #1, SPC(7); USING "Molar Mass Ratio = #.$deet HHH 
kJ/kg"; GI; HM 

PRINT #1, SPC(7); ; "Relative Density = #. eettt #HRHE 
kJ/Ideal CM"; GR: 

PRERT #2). SPOC7) 7 Pactoe coe Sadee aged 
kJ/Real CM"; Z; HVR 

GOTO 340 

320 HM = HM * .4299226# 

HV = HV * .02683919# 

HVR = HV / Z@ 

PRINT SPC(7); USING “Molar Mass Ratio = Fae REE SHERE 
Btu/lb"; GI; HM 

PRINT SPC(7); USING “Relative Density = RLARRES HHH # 


Bru/Ideal CF"; GR; HV 
PRINT SPC(7)}; USING "Compressibilty Factor = #.##### : HHH LH 


Btu/Real CE": 2: HYR 


Btu/Real 7 HYVR 


TF FS = "¥" THEN 330: TF FS = "y" THEN 330 

GoTo 350 

330 PRINT #1, : PRINT #1, SPC(7); USING *Molar Mass = #H. HE 
GROSS HEATING VALUES"; RM 


PRINT #1, SPC(7); USING "Molar Mass Ratio = Fate HHHES 
Btu/lb"; GI; HM 
PRINT #1, SPC(7); USING "Relative Density = RPEeRE HHH # 


Btu/Ideal CF"; GR; HV 
PRINT #1, SPC(7); USING "Compressibilty Factor = #.##### HEHE E 
Btu/Real CF"; 2; HVR 
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340 PRINT #1, 
FESSOR ICT ICI IOGIIUIOUISIGOIOIDOOUISIOIOICISICOIOI IOI UUICOIG ICI ICR ICO I IG IGE! 
350 PRINT 
eee TT TET ETT E TTT TEST TECTCTTTCTT TCT SETTTE TEESE CT CT TTT CCT TTT Tee ere errer rn 
REM 
REM **#*k Axe RHRE EH TNDUTTING POR HARD COPY ** 4 8k eee a eRe eRe REE EERE 
REM 
PRINT : PRINT "DO YOU WANT A HARD COPY NOW (ONLY APPLICABLE FOR DAT MATRIX OR 
PCL PRINTER)?” 
INPUT Y¥1$ 
i IF y1$ “Y¥" THEN 440 


"y" THEN 440 


360 IF ITR = 1 THEN 590 

IF X(15) = 1 THEN 380 

INPUT "Do you want to calculate water saturated values for custody transfer?", 
c 

IF Y$ = "Y¥" THEN 370 

IF ¥$ = "y" THEN 370 

GOTO 510 


370 XW = X(15): XK{15) = 0 
ZL = ((4,4412543D-13 * TL 
2 


- 8.4150417D-10) * T1 + .0000178383014) 
22 = {Zl * Pl - .028354721#) * Ti + 18.87643854# - (2991.2729# ¢ Ti + 
6017.0128#) / T1 
PW = EXP(Z2 + 2.858497# * LOG(T1)} 
FOR I= 1 TO 21: X{I) = X(T} * (1 - PW / PP) / (1 - XW): NEXT I 


X(15) = Pw / PP 

HN = (1 - X(15)) * BNW 

380 ITR = ITR +1 

CLS : PRINT : PRINT ; PRINT : PRINT : PRINT SPC(35); "SATURATED ANALYSIS": PRINT 
IF FS = "Y" THEN 390 

IF FS = "y" THEN 390 

GOTO 280 


400 IP Fg = "¥" THEN ¢ 


s 
fc] 
PRINT ; PRINT : PRINT : PRINT "DO YOU WANT TO DO ANOTHER CALCULATION? " 


INPUT Y$ 
IF ¥$ = " THEN 410 
IF Y$ = “y* THEN 410 
GoTo 490 


410 INPUT "DO YOU WANT TO SAVE OUTPUT TO A FILE"; FS 

IF FS = "Y" THEN 420 

IF FS = "y" THEN 420 

GOTO 430 

420 PRINT "What drive and filename would you like to save the output on?” 
PRINT “For example, B:STREAML.DAT" 

INPUT "What is the desired drive and filename"; FileNames 

OPEN "o", 1, FileName$ 

430 PRINT : PRINT : PRINT "DO YOU WANT TO CHANGE THE TEMPERATURE OR PRESSURE?" 
INPUT K$ 

PRINT : PRINT "DO YOU WANT TO CHANGE YOUR MOLE FRACTIONS!" 

INPUT LS 

IF L§ = "¥" THEN 490 
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IF L$ = *y" THEN 40 

IF K$ = "Y" THEN 40 

IF KS = "“y" THEN 4g 

IF KS = "y" THEN 40 

GOTO 8U 

REM 

REM +#+# eee PRINTING HARD COPY (BUGS OEE SOOO ODE OEE eEe 
REM 


440 IF ITR <> 1 THEN LPRINT SPC(21); “NATURAL GAS PROPERTIES FROM COMPOSITION": 
LPRINT : LPRINT 
IF TA >= 30 THEN 450; IF TB >= 30 THEN 450 


IF ITR <> 1 THEN LPRINT "“CONDITICNS:"; SPC(1); "T(heating value} ="; TB; "DEG 
C*; SPC(2); “Tidensity) =*; TA; "DEG C"; SPC(2}; "P ="; P; "KPA": LPRINT 

GOTO 460 

450 IP ITR <> 1 THEN LPRINT "CONDITIONS: ", SPC(1); "T(heating value) Bs 
"DEG FP"; SPC(2); "T(density) ="; TA; "DEG F"; SPC(2}); "BP ="; PB; "PSIA"; LPRINT 
460 TR TPR <> 1 THEN TPRINT SPC(35); "ANALYSTS": LPRINT — 


IF ITR2 = 1 THEN LPRINT SPC(33); "WET ANALYSIS" 

IF ITR2 = 1 THEN LPRINT "CONDITIONS LPRINT : LPREINT "RELATIVE HUMIDITY ="; 
BA; SPC(3); “HUMIDITY OF GAS 7; AG; SPC(3); “FRACTION OF EXCESS AIR ="; EA 
IF ITR = 1 THEN LPRINT SPC(30); "SATURATED ANALYSIS": LPRINT 

FOR I= 1707 


LPRINT SPC(7): USING "x\ \= ote xX \= #. FoR x\ \= 
H.eHHHG NAS(T); X(T); NAS(I + 7); X(I + 7); NAS(I + 14); X(I + 14) 

NEXT I 

LPRINT : LERINT SPC(7); USING "Molar Mass = HE ai GROSS 


HEATING VALUES"; RM 
IF TA >= 30 THEN 470: IF TB >= 30 THEN 4706 


LPRINT SPC(7); USING “Molar Mass Ratio = $.ehbet THREE 
kJ/kg"; GT; HM 


LPRINT SPC(7); USING "Relative Density = F.REHRE TE 
kd/Ideal CM"; GR; HV 

LPRINT SPC(7}: USING "Compressibilty Factor = #.####¢4 HREEH 
kd/Real CM"; 2; HVR 

470 LPRINT SPC(7); USING "Molar Mass Ratio = F eetee HEEES 
Btu/lb"; GI; HM 

LPRINT SPC(7}; USING "Relative Density = FREE att oF 
Btu/Ideal CF" R; HV 

EPRINT SPC(7}; = Faaatt RHEE LF 


Btu/Real CP"; 2; HVR 

IF ITR2 = 1 THEN GOTO 480 

LPRINT 

BUI ISIS IO IGG ICICI IOI TOG IIICIGEIDIO ORIG SI IOI TE IOIGIGSI GRIGIO a gk dod de dedi 
480 IF ITR = 1 THEN 500: IF ITR2 = 1 THEN 400 

GOTO 360 

490 END 

500 IF ITR2 = 1 THEN 400 

REM 

REM CCGG IR GnG GG ooGdbE RG WED CALCULATIONS (2001 Gtr titi 
REM 

510 INPUT "De you want wet values for engineering calculations?", Y$ 

ITR = 1: [TR2 = 1 

IF Y$ = "Y¥" THEN 520 

IF ¥$ = "y” THEN 520 

GOTO 400 

520 INPUT "What is the relative humidity of air (use fraction not percent)?", HA 
PRINT : PRINT "Enter a 0 for relative humidity of gas to default to X(water) 
inputted earlier.": PRINT 
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INPUT “What is the relative humidity of the gas (use fraction not percent)?", HG 
INPUT "What is the fraction of excess air (use fraction not percent)?", EA 
FOR T = i TO 21 

X(I} = MA(I) 

NEXT I 

ALPHA = 0: BETA = 0: ITR2 = 0 

FOR I = 1 Ta 4 

ALPHA = ALPHA + I * X(I) 

BETA = BETA + (3 + I + [I - 1)) * X(T) 

NEXT 1 

FOR IT = 5 TO 6 

ALPHA = ALPHA + (I - 1) * X{Z) 

BETA = BETA + 2 * I * X{I) 


NEAT I 

FOR I= 7 TO 8 

ALPHA = ALPHA + {I - 2) * K(Z} 
BETA = BETA + (2 * I - 2) * X(I) 
NEXT I 


BETA = BETA + 2 * X{19) 

GAMMA = X{19} 

XW = X{15): K(15) = 0 

Z1 = ({4.4412543D-13 * Tl - 8.4150417D-10) * T1 + .000017838301#) 

22 = (21 * Tl - .028354721#) * TL + 18.87643854# - (2991.2729# / T1 + 


87# * LOG(T1)) 


IF ABS(HGA - HG) > .2 * HG THEN PRINT “INCONSISTENT VALUE FOR HG AND X(WATER)" 
FOR I = 1 TO 21 


X(I) = (1 - X(15}) * X¢Z) 

NEXT I 

HN = (1 - X(15)) * HNA 

NGW = (HG * PW / PP) / ({1 - X(17) - K(21)) * (1 - HG * PW / PP)} 


NAW = 4.77418 * (ALPHA + BETA / 4 + GAMMA) * (1 + BA) * (HA * PW / PP) / (1 - HA 
* PW / PP) 


N = (ALPHA + GAMMA + (X(17) + X(21)) / (1 - X(L7) - K(21}) + (ALPHA + RETA / 4 
+ GAMMA) * (3.77418 + 4.77418 * EA)) 
HN = HN + (NGW + NAW - NuW * (pW / PP) / (1 - PW / PP}) * H(i5) 


ITR2 = ITR2 + 1 

CLS : PRINT : PRINT : PRINT : PRINT : PRINT : PRINT SPC(33}; 
PRINT 

PRINT : PRINT "CONDITIONS:": PRINT 


PRINT "RELATIVE HUMIDITY 7 HA; SPC(3}; "HUMIDITY OF GAS ="; HG; SPC(3}; 
"FRACTION OF EXCESS AIR ="; EA: PRINT 

IF F§ = "Y" THEN $30: IF FS = "y" THEN 530 

GOTO 280 

530 PF fl, : PRINT #1, SPC(33); “WET ANALYSIS" 

PRINT #1, "CONDITIONS:": PRINT #1, 

PRINT #1, "RELATIVE HA; SPC(3); "HUMIDITY OF GAS ="; HG; SPC(3); 


“FRACTION OF EXCESS AIR 
GOTO 280 


: EA: PRINT #1, 
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